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Abstract 

Ql| We calculate the limits on the fraction of viable dark matter minihalos in the early universe 

to host Population III. 1 stars, surviving today as dark matter spikes in our Milky Way halo. 
Motivated by potential hints of light dark matter from the DAMA and CoGeNT direct dark 
matter searches, we consider thermal relic WIMP dark matter with masses of 5, 10, and 20 
GeV, and annihilation to n + n~ , t + t~ , and qq. From this brief study we conclude that, if dark 
matter is light, either the typical black hole size is < lOOM© (i.e. there is no significant Dark 
Star phase), and/or dark matter annihilates primarily to or other final states that result 

in low gamma-ray luminosity, and/or that an extremely small fraction of minihalos in the early 
universe that seem suitable to host the formation of the first stars actually did. 

1 Introduction 

The very first generation of stars, known as Population III. 1 , likely formed from the pristine gas at 
the centers of ~ 1O 6 M dark matter minihalos at z > 10 [1J. The response of the dark matter in 
a minihalo to the formation of a compact baryonic object at its center is a contraction of the dark 
matter density profile in and around the object. When the first stars, expected to be > 1OOM , 
ended their lives by collapsing to black holes3, each remnant remained surrounded by a region of 
enhanced dark matter density, which we call a dark matter spike. In fact, if the first stage of stellar 
evolution is a Dark Star phase, during which the star is powered by dark matter annihilations, 
the first stars would have grown to be even larger, leaving correspondingly larger black holes and 
surrounding dark matter spikes. 

Many of these spikes, remnants of the formation of the first stars, may have fallen into our 
Galactic halo, constituting Milky Way dark matter substructure today. If the dark matter in the 



1 Stars in the mass range ~ 140 — 260Af Q would have ended their lives as pair instability supernovae, leaving no 
remnants [2]. We do not consider these objects here. 
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spikes is made of Weakly Interacting Massive Particles (WIMPs), typically these are their own 
antiparticles, and they annihilate with one another inside the spikes surrounding the black holes. 
This paper focuses on the gamma-ray flux from the annihilation, and compares it to data from 
the Fermi Gamma-Ray Space Telescope (FGST). In previous work 0131, hereafter SDFS, we 
considered WIMPS with masses in the range 100 GeV to 2 TeV and used the FGST data [|5l[6l 
to constrain models of Population III. 1 star formation and/or dark matter annihilation for these 
WIMP masses. Subsequently, the possibility that WIMP dark matter is non-thermal in origin was 
also investigated [7]. In this paper, we return to the thermal WIMP paradigm, but consider WIMP 
masses of 5, 10, and 20 GeV. 

The motivation for this study is the recent increase in attention to lower mass WIMPs due 
to possible hints of their discovery in direct and indirect detection experiments. The DAMA and 
DAMA/LIBRA direct detection experiments [8 J have for ten years found an annual modulation BH 
of their signal compatible with ~10 GeV WIMPs, and now the CoGeNT 001 experiment claims 
a 2.8cr annual modulation signal, expected to be compatible with that observed by DAMA flTTTl . 
Since the two experiments are made of different detector materials (the relevant nucleus for DAMA 
is Na while for CoGeNT it is Ge), and the two experiments are on different continents, the situation 
is interesting. Additionally, the CRESST-II experiment, using calcium tungstate (CaWC^) as the 
target material, has recently announced an excess of events that also seems to be compatible with ~ 
10 GeV WIMPs [12]. Yet, the CDMS experiment, also made of germanium, claims to exclude the 
range of WIMP masses and cross sections indicated by the annual modulation signal 03l . as does 
the XENON 100 experiment 031 . Further potential evidence for light WIMP dark matter has also 
come from indirect dark matter searches in the form of an apparent excess in gamma-rays from the 
Galactic center region 051 . while more recent gamma-ray studies of dwarf galaxies indicate that 
light WIMPs with the canonical thermal annihilation cross section are disfavored 061 . Certainly 
the experimental situation is unresolved, but still it is interesting to consider the possibility that 
dark matter is made of lower mass WIMPs. Here we examine the dark matter annihilation signal 
from light WIMP dark matter in spikes in our Galactic Halo and compare the expected signal to 
gamma-ray data from the Fermi satellite. 

Several other groups have previously studied signatures of annihilation in dark matter overden- 
sities, or spikes, around black holes 07] - l22l (see also the review in (31). Our work is similar in 
spirit and makes the additional step of direct comparison with FGST data. 



Population III. 1 stars likely formed from metal-free, molecular hydrogen-cooled gas at the center 
of dark matter minihalos. We use the parametrization of Ref. [|23l for the minimum halo mass in 
which star formation could occur: 



We take the maximum halo mass for Population III. 1 star formation to be = 1O 7 M , though, 

due to the hierarchical nature of structure formation, our results are not sensitive to this choice. 



2 The First Stars and Their Dark Matter Spikes 
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At some time between the beginning of Population III. 1 star formation and the end of reion- 
ization at z ~ 6, massive Population III. 1 (and Dark Star) formation must have given way to 
subsequent formation of less-massive stars [|24|. however there are few constraints on when this 
transition occurred. Here, motivated by the work in Ref. [|24|. we consider two scenarios for the 
termination of Population III. 1 star formation; at redshifts Zf = 15 and 11. For each case, we as- 
sume that Population III. 1 star formation was possible in any minihalo with a mass between M%£° 
and at redshift z > z f . 

It is unknown how many minihalos meeting the above criteria actually hosted Population III. 1 
stars. We therefore parametrize the fraction that did as /°. The comoving number density of dark 
matter spikes as a function of redshift is then 

N sp ( Z ) = /° (1 - fmergedUl M BH )) N halo (z), (2) 

where Nhaio(z) is the comoving number density of minihalos in which Population III. 1 star forma- 
tion was possible, and f merged, itself a function of /° and the black hole mass Mbh, is the fraction 
of dark matter spikes to have been destroyed in black hole mergers. In SDFS, we estimate that 
f merged is at most 1/2, and is only significant for the largest black holes and /° ~ 1. For small /°, 
f merged -* and the number density of dark matter spikes is determined only by /°. In the fol- 
lowing analysis, we parametrize the fraction of seemingly-capable minihalos to host a Population 
III. 1 star and survive to become part of the Milky Way substructure today as f s . 

Assuming f s = 1, i.e. that each viable minihalo hosted a Population III. 1 star and that /merged 
is negligible, the z = distribution of dark matter spikes throughout the Galactic halo is obtained 
from the Via Lactea II (VL-II) cosmological N-body simulation [25J (see SDFS (3l|4l for details). 
We note that VL-II does not include the role of baryons in formation of the Galaxy (at epochs much 
later than the initial formation of the first stars and their resultant black holes) which could affect 
the numbers of spikes, most likely by contracting the dark matter halos and therefore also the spike 
distribution. However, the extent of this contraction, if it exists at all (e.g. the dark matter halo 
may actually expand during galaxy formation Il26ll ), is uncertain. We also note that dark matter 
spikes may be destroyed when two minihalos (each containing a black hole) merge, forming a 
close binary. While this process is not included in our model, Ref. (3]| includes an estimate of 
this effect: Mergers would change the number of dark matter spikes in our Galaxy today by at 
most a factor of 2 for the most massive black holes considered, and very little for low mass black 
holes. Finally, dark matter subhalos near the Galactic center may experience tidal disruption, but 
the extent to which this affects dark matter spikes surrounding black holes is uncertain |fT9l . In 
Ref. [3], it was found that constraints on f s are generally robust with respect to uncertainties in the 
distribution of spikes near the Galactic center. We return to this point in section HI 

In Fig. Q] we show the number densities of dark matter spikes inside the Milky Way halo as 
functions of Galactic radius for Zf = 15 (red) andzj = 11 (blue). For larger Zf, Population III. 1 
star formation terminates earlier, so there were fewer stars, and therefore the fewer black holes and 
surviving density spikes today; we find 7983 spikes in our Galactic halo for Zf = 15 and 12416 
spikes in our Galactic halo for Zf = 11. In a similar analysis, Ref. [19] found ~ 1027 spikes in 
our Galactic halc@. For comparison, the total dark matter density profile at z = in VL-II is also 

2 As discussed in SDFS, our approach for finding relevant dark matter minihalos can be contrasted with that of 
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shown (black); although the normalization of these points is arbitrary, it is useful to illustrate that 
the total dark matter profile is more extended than the distribution of black holes with dark matter 
spikes. 

The large population of dark matter spikes in our galaxy is consistent with the existence of 
high redshift quasars. The Sloan Digital Sky Survey has identified high redshift quasars at z ~ 6, 
each presumably powered by a super massive black hole (SMBH), with a mass on the order of 
10 9 M Q . These SMBHs grew from seed black holes [|27l : In the standard picture, the seed black 
holes formed from standard Population III stars, dark stars, or direct collapse models (See [28] and 
references therein). 

Simulations indicate that most seed black holes grew by accreting baryons (see [29] and 11301 
for recent simulations), a scenario consistent with the formation of dark matter spikes. In fact, most 
seed black holes failed to become SMBH and end up as intermediate mass black holes (IMBH), 
with masses in the range 10 3 M to 10 5 M . Simulations also indicate that a large fraction of the 
seed black holes only grew by accretion and never merged. Hence, dark matter spikes may well be 
ubiquitous. 
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Figure 1 : The number density of black hole spikes in the Milky Way as a function of Galactic 
radius for star formation models with Zf = 15 (red) and Zf — 11 (blue) as described in the text. 
We assume f s = l. The black points illustrate the total dark matter density profile at z = 0. 

The density profile of an individual dark matter spike surrounding a black hole is a key issue 
in computing the annihilation signals. In this paper we make the best simple estimates we can 

Ref. lH9l . At z = 18, they populated halos that constituted 3cr peaks in the smoothed primordial density field with 
seed black holes of initial mass 100 M Q . Using an analytical model of halo evolution, they simulated 200 statistical 
realizations of the growth of a Milky Way-sized halo. Instead, we use one Galaxy-mass halo from a very high resolu- 
tion cosmological simulation ( VL-II) to follow potential Population EH. 1 star-forming minihalos to z = 0. See (3) for 
further details. 
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and secondly demonstrate what we consider to be the worst case scenario. The progenitors of the 
black holes discussed in this paper formed initially in the following way: as a protostellar cloud 
started to collapse at the center of the dark matter halo, dark matter was gravitationally pulled into 
the collapsing object along with the baryons. The initial density profile at the center of the halo 
then steepened due to this growing gravitational potential. We treat this gravitational dark matter 
enhancement with adiabatic contraction of the dark matter halo around the central mass. A very 
simple prescription for the adiabatic contraction is the Blumenthal et al. method OTA , which only 
considers halo particles on circular orbits; equivalently, only their angular momentum is conserved 
as the halo is compressed. We use this prescription in our simple estimates. However, early 
work on dark stars met with some skepticism because of exactly this issue. Consequently several 
groups have verified that this simple Blumenthal et al. method is accurate to within a factor of 
two for the case of dark stars. Iocco et al. Il32l (see Figure 1 in their paper) used the Gnedin et 
al. method lT33l and found basic agreement with the enhanced dark matter profiles obtained for 
these early stars using the Blumenthal et al. method. Similarly, Natarajan, O'Shea, and Tan |T3"4l 
also used the Gnedin et al. method and agreed. To address these concerns yet more carefully for 
the case of early stars forming at the centers of minihalos, Ref. Il35l performed an exact calculation 
using the Young method [|36l . which takes into account the conservation of the radial action as 
well as angular momentum (the adiabatic invariants) and again reproduced the basic results found 
using the Blumenthal et al. method (in some cases finding densities higher by a factor of two 
and in others lower by a factor of two; see Figures 3 and 4 in their paper. jE Since there are 
many uncertainties and parameter choices here (e.g. the concentration parameter, the redshift of 
the object at the time of its formation, the WIMP mass (which determines the baryon density at the 
time of dark star formation), the cutoff radius inside which we should assume the dark matter has 
already annihilated away by the time the black hole forms (see below), the size it grows to, etc.), 
we cannot be more precise than this factor of two. Hence we present two sets of results: First, we 
use the Blumenthal et al. method to obtain constraints. Second, we present a worst case scenario, 
namely, we obtain results under the assumption that spike densities are lower by a factor of two 
relative to the Blumenthal et al. prediction. The latter analysis results in lower annihilation signals 
and far weaker bounds, and in this sense these are the most conservative results. 

A second major source of uncertainty in the dark matter profile is the choice of initial dark 
matter density profile prior to adiabatic contraction. As our canonical case, we begin with Navarro, 
Frenk, and White (NFW) profiles ll38l for both the baryons and dark matter, and, as described in the 
previous paragraph, use the Blumenthal et al. prescription for adiabatic contraction OTTl . However, 
it is possible that such an NFW profile overestimates the amount of dark matter at the core of a dark 
matter minihalo where a Population III. 1 star formfl To address this concern, Ref. 11331 looked at 
the extreme case of an initial halo profile that has a core of constant dark matter density, known 
as a Burkert profile [39]. We note that this profile is not realistic for early minihalos and invoke 

3 It is important to point out that Youngs adiabatic prescription is not at all sensitive to departures from sphericity: 
Sellwood and McGaugh [37] showed that, for what concerns adiabatic contraction, even a fiat disk could be well 
approximated as a sphere. 

4 We do note, however, that it is not unreasonable to imagine that in these pristine early halos the central profile 
was steep while subsequent baryonic processes lowered the dark matter central cusp by the present epoch where 
observations are made. 
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it only as a scenario that results in the absolute minimum core dark matter density. Even in this 
extreme case, the dark matter density at the core of the star forming in the center of the minihalo is 
only lower by about a factor of two. We note that an initial Einasto profile H51 (in lieu of NFW or 
cored) would result in a spike density in between these two extremes. Thus we feel confident that 
we may use the NFW case as a good estimate and present results as well for the worst case where 
the dark matter density is lower by a factor of two due to the different possible initial profiles. 

In short, to obtain the dark matter density profile around the black hole we use the Blumenthal 
et al. method to obtain simple estimates of adiabatic contraction starting from an initial NFW halo. 
We then also present results with densities lower by a factor of four: a factor of two due to the 
initial dark matter profile of the minihalo and another factor of two due to the success of adiabatic 
contraction. This factor of four represents our sense of the most extreme possible uncertainties in 
our results. The sensitivity of our results to these choices is discussed in section HI as well. 

Finally, after adiabatic contraction, we impose a high density cut-off at small radii in the spike 
to account for the annihilation of the dark matter throughout the lifetime of the spike; 



where m x is the mass of the dark matter particle, (av) its annihilation cross section times velocity, 
and t BH is the lifetime of the central mass, roughly 1 .3 x 10 10 years for a star that formed atz= 15. 

While standard Population III. 1 stars are expected to have masses of ~ 1OOM , if dark matter is 
capable of self-annihilating, then the first stars may have been powered by dark matter annihilations 
for some period of time prior to nuclear fusion. This first phase of stellar evolution is known as the 
Dark Star phase [41 J and may have lasted anywhere from a few hundred thousand years to millions, 
or even billions, of years. During the Dark Star phase, the star remains cool enough to continue 
to accrete baryonic matter, and may grow to ~ 1OOOM 11421 or even as large as > 10 5 M Q [|43l , 
depending on the details of how the dark matter is depleted and replenished in the star. We note 
that the WIMP mass does affect the final mass of the star, as the dark matter heating is inversely 
proportional to the WIMP mass, however the difference in the final stellar mass is less than a factor 
of two for 100 GeV and 1 GeV WIMPs assuming they move on circular orbits and dark matter is 
depleted in a simplistic way [44]. If the dark matter particle orbits are more complicated, lighter 
WIMPs would result only in a cooler and more extended star during the Dark Star phase [|43l . 

When the dark matter fuel inside a Dark Star runs out, it collapses and heats up to become 
a standard fusion-powered star, which, at the end of its life, will likely undergo core collapse 
leaving a black hole remnant. We use the potential existence of the Dark Star phase, as well as 
the possibility of direct collapse of very massive gas clouds to black holes |[28ll4"5ll , to motivate 
consideration of black holes as large as 1O 5 M . 

3 Gamma Ray Signal from Dark Matter Annihilations 

For a Majorana dark matter particle with mass m x and annihilation cross section times velocity 
(crv), the rate of WIMP annihilations in a dark matter spike is 



m 



Pmax 



{av)t B H 



(3) 




(4) 
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where pdm{t) is the dark matter density as a function of spike radius, with r min and r max defining 
the volume of the dark matter spike in which annihilations occur. 

We choose as a benchmark scenario (av) = 3 x 1CT 26 cm 3 s _1 , in agreement with the measured 
dark matter abundance today for thermal WIMP dark matter, and consider several WIMP candi- 
dates, defined by mass and annihilation channel. We consider annihilations of WIMPs with masses 
of 5, 10, and 20 GeV to Standard Model final states r+r~, and dd. Note that these WIMPs 

are not heavy enough to annihilate to W + W~. While all but the lightest WIMPs considered here 
are heavy enough to annihilate to bb, the resulting photon spectra from annihilation to all quark 
pairs (qq) are nearly identical. We compute the spectrum of photons from annihilation to final state 
/, dNf/dE, with PYTHIA ll46ll . except for the final state n + n~ , in which case the photon spectrum 
comes only from final state radiation [47J. In the following analysis, we assume that the branching 
fraction to each final state is 1, though, in principle, some combination of final states is possible. 

The differential flux of neutral particles from annihilations to final state / in a dark matter spike 
with radius r max located some distance D from our Solar System is given by 

d$ f T dN f 

— — = — (5) 

dE AtiD 2 dE" W 

for D 3> r max . If D is not much larger than r rnax , and the spike is not point-like, an integral along 
the line-of-sight to the spike must be performed. 

If a single spike is a sufficiently bright and compact source of gamma-rays, it may have been 
identified as a point source and recorded in the FGST First Source Catalog [5], or, indeed, in the 
EGRET Source Catalog [48] . Spikes bright enough to be identified as point sources must not be 
brighter than the brightest source in the the FGST catalog. Requiring that the flux not exceed 
this brightness establishes a minimal distance, D^f n , beyond which the spike must be located. 
Similarly, we can define a maximal distance, D^ ax , as the distance beyond which a spike would 
not be bright enough to have been detected at 5cx significance (see and [|49ll for details). All 
spikes which could have been identified as point sources lie between D^f n and D^ ax . 

In the left panel of Fig. [2} we display the minimal distance, D^f n , at which a point source may 
be located in order not to exceed the largest flux from any point source measured by FGST for 10 
GeV WIMPs annihilating to dd (green dotted), t + t~ (blue dashed), and (magenta solid). 

For 10 GeV WIMPs annihilating to dd or t + t~, if the central black hole is < 100Af o (roughly 
what is expected from standard Population III. 1 star formation), spikes may be located within 1 
kpc of our Solar System. For annihilations to spikes may be < 1 kpc from our Solar System 

even for black holes as large as C(10 3 )M Q . However, for either choice of Zf considered here, the 
VL-2 simulation results indicate less than one dark matter spike within 1 kpc of our Solar System, 
even for f s — 1. Finally, we remind the reader that we consider only potential FGST point sources 
in this analysis; if dark matter is very light (m x < 5 GeV) and annihilates to yU + /i~ in the spike 
surrounding a very small black hole (M BH < 1OM ), if the nearest spike is actually as near as 
D^f n it may appear to be an extended gamma-ray source for FGST. We do not consider extended 
sources here, though we return to this point in section HI 

The right panel of Fig. [2] shows the maximal distance, D^ ax , at which a single dark matter 
spike would have been identified by FGST as a > 5cr point source. Sources that are farther from 
our Solar System than D^ ax would be too faint to have been identified as point sources in the 
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Figure 2: In the left panel, we display the minimal distance from our Solar System of a single spike 
such that it does not exceed the flux of the brightest source in the FGST First Source Catalog, D^f n , 
as a function of central black hole mass for Intermediate Zf. From top to bottom, the contours are 
for 10 GeV WIMPs annihilating to dd (green dotted), t + t~ (blue dashed), and (magenta 
solid). In the right panel, we show the maximal distance from our Solar System of a single spike 
such that it would have appeared as a > 5a point source to Fermi in the first year of data, D^ x , 
for the same cases. The horizontal black dashed line indicates the distance to the Galactic center 
from our Solar System. 



first year of FGST operations, and therefore the gamma-rays from these spikes would be part of 
the measured diffuse gamma-ray flux. In fact, for 10 GeV WIMPs annihilating to if the 

typical black hole size is quite small (~ 1OM ), spikes within only a few kiloparsecs of our Solar 
System would be too faint to have been identified as gamma-ray point sources by FGST In all 
cases discussed here, for all but the largest black holes considered, at least some spikes in the 
Milky Way halo would contribute to the diffuse gamma-ray flux. 



4 Constraining f s 

As alluded to in section [2l it is quite possible that / s <C 1. In this section we discuss how gamma- 
ray data can be used to constrain f s , given a star formation scenario and a dark matter model. 

First, it is possible to constrain f s with the measured diffuse gamma-ray flux: Since the diffuse 
gamma-ray background is not yet well-understood, we conservatively require only that the total 
diffuse gamma-ray flux from dark matter annihilations around spikes in the Milky Way halo not 
exceed the measured flux in any of the nine FGST energy bins in Ref. [6] by more than 3a. In 
SDFS, the diffuse gamma-ray flux measurements of FGST were found to constrain f s mainly for 
dark matter annihilating to leptonic final states for WIMP masses > 100 GeV. However, for the 
lighter WIMP masses < 20 GeV considered in this paper, no interesting bounds result. Because 
the FGST measurements of diffuse flux indicate that it scales roughly as E~ 2A , diffuse gamma-ray 
constraints on dark matter annihilation tend to come from the highest energy bins, if they exist at 
all. For light WIMPs (relative to standard 100 GeV WIMPs), the entire gamma-ray spectrum is 
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shifted to lower energies such that the gamma-ray flux predicted for WIMPs annihilating in dark 
matter spikes is not in conflict with FGST measurements. Hence we restrict our discussion in the 
remainder of the paper to point sources observed by FGST. 

A second way to constrain f s , making use of FGST point source data, is to require that the 
integrated luminosity of the brightest (i.e. closest) dark matter spike not exceed the integrated 
luminosity of the brightest Fermi point source. As discussed above, the left panel of Fig.|2]displays 
this minimal distance, D^ n , for 10 GeV WIMPs. A limit on f s may be established by requiring 
the expectation of finding less than one spike within a sphere of radius r = D^ n centered on our 
Solar System. That is, 



where N sp (R, f s ) is the number density of spikes as a function of Galactic radius, R, and the 
fraction of spikes to survive, f s . As discussed in section [2l we assume that the actual number 
density of spikes, N sp (R, f s ), is related to the total possible number density of spikes, N sp (R, f s = 
1), as 



In Figure |3l we plot the maximal allowed value of f s as a function of black hole mass for 
Zf = 15 (red circles) and Zf — 11 (blue squares). The filled points correspond to the constraint 
derived from the requirement that the brightest spike not be brighter than the brightest FGST source 
(corresponding to D%? n as plotted in Fig.|2]), which is known to be associated with the Vela pulsar. 
Although we cannot exclude the possibility that the brightest dark matter spike is located along our 
line of sight to Vela, it would be a bizarre and grand conspiracy if all of the brightest dark matter 
spikes were to be located along our lines of sight to the brightest gamma-ray point sources, many of 
which are associated with known astrophysical objects. Thus, alternatively, one could consider the 
constraints derived from the requirement that the brightest dark matter spike not be brighter than 
the brightest unassociated FGST point source; that is, a point source that is not associated with any 
known astrophysical object. These less conservative, but perhaps more reasonable, constraints are 
shown as the open points in Figure [3] 

We also show, as solid green triangles in each panel of Figure[3l the constraints resulting from a 
dark matter density spike that is less dense by a factor of four as discussed above Equation[3]as the 
most extreme possibility; for these points, we required that the brightest spike not be brighter than 
the brightest FGST source and show only the star formation history leading to the weakest con- 
straints. They therefore represent the most conservative limits on f s , accounting for uncertainty in 
the contracted halo profile. These limits (solid green triangles) represent the "worst case scenario" 
for indirect detection (among those considered here). We also note that the difference between 
the solid and open blue and red points is a factor of 4.7 in spike distance or, equivalently, spike 
density, so that the green triangles for the case of brightest unassociated FGST point source would 
lie slightly below the solid blue points (again, assuming the less favorable scenario for the termina- 
tion of Population III. 1 star formation). We remind the reader that throughout we have assumed a 
thermal annihilation cross section, (av) = 3 x 10~ 26 cm 3 s -1 ; a larger (smaller) annihilation cross 
section would lead to stronger (weaker) constraints than those shown in Figure [3] 

In all cases, constraints are strongest for the largest black holes, which result in the most dense 




(6) 



N sp (R, f s ) = f s x N sp (R, f s = l). 



(7) 
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Figure 3: Upper limits on the fraction f s of minihalos to host the first stars, plotted as a function of 
the mass of the remnant black hole at the center of each dark matter spike. These limits are derived 
from the gamma ray flux from point sources measured by FGST. Squares and circles correspond 
to our best estimates while the green triangles illustrate the extreme case of most conservative 
bounds as described in the text. The filled points are the upper limits obtained by requiring the 
dark matter annihilation signal in gamma-rays not exceed the brightness of Vela, while the open 
points illustrate the upper limits obtained by requiring the dark matter annihilation gamma-ray 
flux not exceed that from the brightest unassociated FGST point source. Intermediate and Late Zf 
are shown as red circles and blue squares, respectively. 

spikes, and for the lightest dark matter, for which the number density in a spike is the largest. 
For XX ~^ l^l 1 ' ■> shown in the top panels, the gamma-rays come exclusively from final state 
radiation, resulting in relatively weak constraints. For example, if the typical black hole mass is 
(9(10 — 100) M , the fraction f s is unconstrained for any dark matter mass > 5 GeV. However, even 
for such light dark matter, the constraints can be quite strong if the typical black hole mass is larger. 
In the middle and bottom panels, for \X ~^ t + t~ and XX ~^ dd, respectively, annihilations lead to 
states that either hadronize or decay to states that hadronize, resulting in a larger gamma-ray flux 
and therefore stronger constraints on f s . Remarkably, in these cases constraints are at the percent 
level for black holes of mass ~ 1OOM , the typical mass expected from standard Population III. 1 
star formation (in the absence of a Dark Star phase). 

As discussed in section [3j we consider only FGST point sources in this analysis. If the dark 
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matter is very light (m x < 5 GeV) and annihilates to in the spike surrounding a very small 

black hole (M BH < 1OM ), and if the nearest spike is actually as near as D^ n , it may appear to 
be an extended gamma-ray source for FGST. In the upper left panel of Fig. |3j it is apparent that any 
scenario in which there is the possibility that the nearest spike might appear as an extended source 
is entirely unconstrained anyway. An analysis of extended sources may prove very interesting: it 
may be possible to discover the distinct shape of a dark matter spike in this manner, though no 
additional constraints would result for the scenarios considered here. 

In principle it might be prudent to exclude the inner few kpc from the Galactic center from our 
study because spikes in the inner region of the Galaxy may have been disrupted, though in practice 
the results do not differ significantly from what we report here. The origin of the insensitivity to 
the Galactic center region can be understood by considering the following: First, models in which 
^min ~ 3.5 kpc are not affected at all by what happens within ~ 5 kpc of the Galactic center. 
Second, for models in which D^ n is comparable to the radius of the Solar Circle, we demand 
less than one spike within a volume of ~ 10 3 kpc 3 centered on our Solar System, already orders 
of magnitude below the expected number density of spikes in the Solar neighborhood, so these 
cases are strongly constrained whether or not the spikes nearest the Galactic center are considered. 
Furthermore, the vast majority of spikes are not located near the Galactic center; for Zf = 15 and 
11, the percentages of Milky Way spikes within 5 kpc of the Galactic center are just 10% and 4%, 
respectively. Thus, the constraints on f s are robust with respect to uncertainties in the dynamics 
near the Galactic center. 

Finally, we would like to stress that it is possible that one or more point sources in the FGST 
catalog are, in fact, dark matter spikes. In this case, spectral information may provide some clues 
as to the nature of dark matter. While a spectral analysis is beyond the scope of this study, Ref. [49 1 
finds that as many as 20 to 60 gamma-ray point sources may in fact be dark matter substructure. 
Perhaps we have more information about dark matter spikes than we imagine. 

5 Discussion and Conclusion 

We have calculated the limits on the fraction of viable dark matter minihalos in the early universe to 
host Population III. 1 stars, surviving today as dark matter spikes in our Milky Way halo. Motivated 
by potential hints of light dark matter from the DAM A and CoGeNT direct dark matter searches, 
we considered thermal relic WIMP dark matter with masses of 5, 10, and 20 GeV, and annihilation 
to t + t~, and qq. We find that constraints from the integrated luminosity of Vela are only 

significant if the typical black hole mass is > 10 3 M Q , as would be expected if the first stars 
experience a Dark Star phase of stellar evolution, during which the star is powered by dark matter 
annihilations. Excluding the very unlikely coincidence that all of the brightest dark matter spikes 
are located along our lines of sight to the brightest known astrophysical gamma-ray sources, then 
the constraints on f s are at the percent level even for 100M o black holes if dark matter annihilates 
primarily to quark pairs or tauons. Even in our "worst case scenario" for indirect detection of dark 
matter spikes, constraints are significant for typical black hole masses > 1O 3 M . 

From this brief study we conclude that, if dark matter is light, either the typical black hole size 
is < IOOMq (i.e. there is no significant Dark Star phase), and/or dark matter annihilates primarily 
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to or other final states that result in low gamma-ray luminosity, and/or that an extremely 

small fraction of minihalos in the early universe that seem to be suitable to host the formation of 
the first stars actually did. We look forward to additional evidence of early star formation and the 
properties of dark matter, and, eventually, to a solid understanding of both. 
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